Saturation of the production of quantum entanglement between 
weakly coupled mapping systems in strongly chaotic region 

Atushi Tanaka[| 

Department of Physics, Tokyo Metropolitan University, 
Minami- Osama, Hachioji, Tokyo 192-0397, Japan 

Hiroshi Fujisaki[] 

Department of Theoretical Studies, Institute for Molecular Science, 
Myodaiji, Okazaki 444-8585, Japan 

Takayuki MiyaderaQ 

Department of Information Sciences, 
Tokyo University of Science, Noda City, Chiba 278-8510, Japan 
(Dated: February 1, 2008) 

Abstract 

The production of quantum entanglement between weakly coupled mapping systems, whose 
classical counterparts are both strongly chaotic, is investigated. In the weak coupling regime, it 
is shown that time correlation functions of the unperturbed systems determine the entanglement 
production. In particular, we elucidate that the increment of nonlinear parameter of coupled kicked 
tops does not accelerate the entanglement production in the strongly chaotic region. An approach 
to the dynamical inhibition of entanglement is suggested. 
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In a quantum composite system, even if the subsystems are remotely separated and the 
whole system is in a pure state, the subsystems generically have a nonclassical correlation [|1|. 
This striking phenomena is called quantum entanglement [Q, which is utilized not only to 
achieve the procedures that have no classical analogs (e.g. quantum information process- 
ing [Q), but also to realize the "classical world" in which quantum interference phenomena 
are "decohered" as a result of quantum dynamics Even when there is no quantum entan- 
glement between subsystems, a weak interaction between the subsystems generally produces 
the quantum entanglement during unitary time evolutions This is an important dynam- 
ical origin of decoherence @]. 

Through a number of numerical experiments, it is known that the productions of en- 
tanglements induced by quantum dynamics heavily depend on the qualitative nature of the 
corresponding classical dynamics, namely, regular or chaotic @, ^, ||, || , as is easily expected 
from the studies of "quantum chaos" |10[- On one hand, in classically regular systems, the 
confinement of phase-space dynamics in a narrow region enclosed by KAM tori make it 



difficult to produce strong entanglements in general ||rT|. On the other hand, the absence 
of such dynamical barriers in classically chaotic systems promotes the production of entan- 
glement. Although there are quantum effects on the phase-space dynamics, e.g. tunnelings 
and localizations |T^, in both regular and chaotic systems, it is confirmed that the scenario 
above quahtatively holds [§|, 0, §], 3- 

This motivates the next question: In the chaotic region, does stronger chaos enhance 
the production of entanglement? Looking for an analogy of a study on quantum open sys- 
tems fl^. Miller and Sarkar obtained a numerical result that suggests the linear instability 
of classical dynamics enhances the production of entanglement 0. Their numerical experi- 
ment however concerns only in the weakly chaotic region where chaotic seas and tori coexist. 
The complexity of phase-space dynamics is the source of the difficulty to obtain a theoretical 
explanation for the Miller and Sarkar 's result. 

Our aim is to provide a theoretical argument of entanglement production in weakly 
coupled chaotic systems. In contrast to the Miller and Sarkar's work, we focus on the strongly 
chaotic region where the effect of tori is small, to facilitate to obtain a theoretical explanation. 
Starting from separable pure states, we examine the productions of quantum entanglement 
due to unitary time evolutions. The entanglement production processes are slow due to the 
weak coupling. Furthermore, the recurrence time of classically chaotic systems is relatively 
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long. Hence, the entanglement production processes are nearly stationary processes, at 
least, in a short time period. This enables us to introduce an entanglement production rate. 
We investigate how the entanglement production rate depends on the nonlinear parameter 
below. 

Our numerical experiments employ coupled kicked tops Firstly, we intro- 

duce their constituent, a kicked top |T3, which is described by the Hamiltonian 



Hk = IT jy / 2 -\- A{t) kJ'^/{2j), where Jj is the i-th component of the angular momentum 
operator of the top, j is the magnitude of the angular momentum, A; is a nonlinear param- 
eter, and A(t) = J2nez^i^ — n) is a "periodic delta function". Secondly, we employ the 
following Hamiltonian to describe the whole system that is composed of two kicked tops: 

H = Hk,®i + l(»Hk, + e VA{t) (1) 

where e is a coupling constant, V = JizJ2z/j is the interaction Hamiltonian, and Jiz is the 
Jz of the i-th top. We report the case that the magnitudes of the angular momenta of the 
two subsystems are the same value, j. 

Since we focus on the case where the total system is in a pure state, our choice of a measure 
of quantum entanglement between the two subsystems, is the von Neumann entropy S^n of 
the first subsystem |T^. Note that the von Neumann entropy of the second subsystems is 



equal to the one of the first subsystem, when the total system is in a pure state. At the same 
time, we employ the linear entropy Sn^ instead of ^vn, to facilitate theoretical arguments. 
In our numerical experiments, these entropies behave qualitatively similarly. 

We examine the productions of quantum entanglements, by observing and Sun, during 
unitary time evolutions whose initial states are product (i.e. separable) states. A typical 
result in the strongly chaotic region is shown in Fig. 0: When the coupling constant e 
is small, there is a t-linear entanglement production region, which is wide enough to call 
it "stationary" entanglement production. Note that during the stationary entanglement 
production, the state vectors of the subsystems are spread over the phase space of the 
subsystems (Fig. ^. In other words, the phase space distribution of each subsystem is nearly 
in "equilibrium". We accordingly expect that the each subsystem plays a role of a chaotic 
"heat bath" of its partner [^. The t-linear entanglement production region starts at a time 
step T', after a short transient to attain the equilibrium of the phase space distribution of 
the subsystems, and ends at a time step T", until the increment of the entropy reaches its 
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FIG. 1: Time evolutions of quantum entanglement, measured by the entropies of a subsystem. 
The entropies are scaled by 5*0 = 2e^j^ (cf. eq. @). From bottom to top, we depict with 
e = 10~^, 5 X 10~^, 3 X 10~^, 10~^, 5 x 10~^, 10~^. An estimation Sy^^ by our perturbation 
formula is degenerate with the case e = 10~^. In the inset, the solid line and the dotted 
line correspond to S\\^ and 5vN) respectively, at e = 10~^. The value of nonlinear parameters 
are fci = /c2 = 7.0, which means that the corresponding classical tops are strongly chaotic |15]. 



The magnitude of the angular momenta is chosen to be large j = 80, in order to investigate the 
semiclassical regime. The center of the initial wave packet, which is a product of spin-coherent 



states |17|, is ((91, 0i, 6*2, (^2) = (0.89,0.63,0.89,0.63). 



equilibrium (see Fig. larger e). 

In order to explain the t-linear, stationary entanglement production, we employ a time- 
dependent perturbation theory, whose small parameter is a coupling constant e, to evaluate 
the linear entropy S'iin(t) at t-th step. The resultant formula is 



t t 



C^w = ^oEE^K^) (2) 



m=l n=l 

.2-2 



where Sq = 2e j , and D(m,n) is a time-correlation function of the uncoupled system. Fur- 
thermore, since the interaction Hamiltonian V takes a bilinear form, D{m, n) is decomposed 
as follows 

D{m,n) = Ci{m,n) C2{m,n) (3) 
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FIG. 2: The Husimi function [18| of the first subsystem at t = 15, during the stationary entan- 



glement production region. The contour and density plots are in normal and logarithmic scales, 



respectively. Note that this region is beyond the Ehrenfest time [19|. We choose e = 10 ^. Other 
parameters are the same as in Fig. |l[ 

where Cj(m, n) = j^"^ {{Ji^ Jiz) — 01^)011)) is a normalized correlation functions of J^" , which 
is evolved by the unperturbed Hamiltonian H\^=o until n-th steps, with an initial condition 
j£. = Jiz, and the expectation value (■) is respect to the unperturbed system. The details 
to obtain the formula (0) will be shown elsewhere 

We remark on the entanglement production formula (0): (i) Although we start form 
the evaluation of the entropy of a subsystem, the formula is in a symmetric form with 
respect to the exchange of the two subsystems. This is consistent with the symmetric 
nature of quantum entanglement when the whole system is in a pure state; (ii) Since our 
approach does not take into account the effect of the recurrence, the formula would have 
qualitatively different applicability to the classically regular and chaotic systems. On one 
hand, for classically regular systems, our theory would break down in relatively short time 
period, due to the smallness of the period of the recurrence. On the other hand, for chaotic 
systems, we numerically confirmed that our theory works rather long time period; (iii) Our 
formula has a similarity with the ones in phenomenological descriptions of linear irreversible 



5 




10 20 30 40 50 
T 



10 



10^ 



+ 



q 10" 



10" 







(b) 



+ a 



B X ^ 

+ □ ± 



10 

T 



15 20 



FIG. 3: The r dependence of the correlation function \D{t + r, t)| for (a) a regular system (ki = 
k2 = 1.0) and (b) a chaotic system (ki = k2 = 7.0) with j = 80. Different symbols correspond 
to different values of t (+: t = 40, x: t = 50, *: t = 60, □: t = 70). Note that (b) employs a 
normal-log scale. 



processes p2[, in the sense that these theories use time correlation functions to describe 
relaxation phenomena. This is useful both for discussing phenomenological arguments and 
for making a link with a phenomenological theory and a microscopic theory p2 |. 

Before applying the formula (§) numerically, let us confirm that the time correlation func- 
tion D{m,n), which is the most important ingredient of the formula (0), strongly depends 
on the dynamics of the classical counterparts (Fig. |^). On one hand, in the regular case, 
D{m,n) decays slowly with large oscillations, as the time interval |m — n| become large. On 
the other hand, the chaotic dynamics makes the decay of D(m, n) much faster. Such a rapid 
convergence of the correlation function, together with the formula (0), implies the t-linear, 
stationary entanglement production region (see Fig. |T]). 

The perturbation formula (0) provides an approximate estimation of the entanglement 
production rate F of the t-linear, stationary entanglement production region: 

J^// rp// 

r^"-^o^;r^E E^K-) (4) 

m=T' n=T' 

where T' and T" are the start and the end of the t-linear entanglement production region, 
respectively. In the strongly chaotic region, where the effect of tori is small, it is possible to 
give an analytical estimation for F^""" (H). Since D{m, n) is the product of the fluctuations of 



Jizi whose distribution functions become quickly uniform due to the chaotic dynamics (see 
Fig. we assume that D{m,n) decays exponentially 

D{m,n) = DQexp{—'y\m — n\) (5) 

The prefactor Dq is determined by the magnitude of the fluctuations of Ji^ and J2z, whose 
distribution functions are almost uniform in strongly chaotic systems (see Fig. §). We 
accordingly assume Dq = (1/3)^, which is independent of ki and ^2. Furthermore, it is 
natural to expect that the decay rate 7 of D{m,n) increase as the degree of chaos of the 
classical counterpart becomes stronger (i.e., as the values of nonlinear parameters ki and k2 
increase), when the effect of tori is negligibly small, although we could not obtain a precise 
nonlinear parameter dependence of 7 due to the fact that the decay is very fast (within a 
step) in the strongly chaotic region (see Fig. 0(b)). 

As long as the t-linear region is wide enough, i.e. T" — T' ^ 7^^, the exponential decay 
assumption (||) provides an estimation 

rPT~5oDocoth(7/2). (6) 

This provides our main result: in the limit that the classical counterpart is the strongly 
chaotic, i.e. in the limit 7 00, the perturbation theory (^ predicts that the entanglement 
production rate F converges to a finite value. Furthermore, the convergence is expected to 
be fast, when 7 is larger than unity. Our main claim is that the effect of the enhancement of 
entanglement due to chaotic dynamics saturates in the strongly chaotic region, in contrast 
with the weakly chaotic systems |0, |^, This prediction is confirmed by our numerical 
experiments (Fig. We explain the reason why the strongly chaotic systems have smaller 
entanglement production rates, in comparison with the weakly chaotic systems. As the 
chaos becomes stronger, the correlation time scale of the interaction Hamiltonian V, which 
is evolved by the unperturbed Hamiltonian Hq in the interaction picture, becomes shorter. 
Hence, due to the dynamical averaging, the infiuence of V is effectively reduced. Conse- 
quently, the entanglement production rate is also reduced. In particular, when the chaos 
is strong enough, the effect of the perturbation on F^""" (P) comes only from the "diagonal" 
part D{n,n). This is the origin of the saturation at large k (Fig. ^). 

We make a brief remark on our result, in correspondence with the existing publica- 
tions that suggests the entanglement production rate is proportional to the Lyapunov ex- 
ponent of the classical counterpart (i.e., chaos promotes quantum entanglement) |14|. 
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FIG. 4: Dependences of the entanglement production rates F, which is measured by hnear entropy, 
on the nonhnear parameter k = ki = k2- In order to show typical k dependences, we choose 
several initial conditions (depicted by different marks) that place in the chaotic sea. Although 
the entanglement production heavily depends on the initial condition in the weakly chaotic region, 
the disappearance of tori weakens the initial condition dependence in the strongly chaotic region. 
Other parameters are the same as in Fig. ^. 

Firstly, although our model system in the numerical experiments is the same as Miller and 
Sarkar's workQ, the result is qualitatively different. The difference comes from the differ- 
ent "strength" of chaos. In contrast with the strongly chaotic systems, the entanglement 
productions of the weakly chaotic systems are significantly influenced by the existence of 
tori. A crudest explanation is provided by the perturbation formula (^: During the chaos 
is not fully developed, a large portion of the phase space is occupied by tori. This reduces 
Dq, which is the magnitude of the fluctuation of the interaction V. Accordingly, as the 
chaos become stronger in weakly chaotic region, the development of chaos increases Dq. 
This results in the increment of the entanglement production rate. In the strongly chaotic 
region, the growth of Dq saturates due to the breakdown of tori and 7 takes a major role in 
the parametric variation of F^^ (^). Secondly, we compare our result with Zurek and Paz's 



work on open systems [Ij], since the each subsystem in our numerical experiment acts as 
"heat bath" of its partner. The most important difference comes from the fact that the two 
works focus on completely different region, far before (Zurek and Paz) and far after (ours) 
the Ehrenfest time (see. Fig. H). 

We expect that our result on the saturation of entanglement between weakly coupled 
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systems is generic, in strongly chaotic systems with a "compact" phase space that allows 
the assumption (H). However, we note that our result heavily depends on the fact that 
the system is a discrete-time system, i.e. a mapping system. 

Finally, we discuss an extension of our result to flow systems. It is straightforward to 
obtain a perturbative estimation of entanglement production rate F oc I/7, which implies 
the suppression of the entanglement productions in the strongly chaotic limit 7 — 00. This 
is completely opposite to the case in the weakly chaotic region |P, 0, |, If. More thorough 
investigations on this point will be reported in subsequent publications. We remark that a 



similar suppression of quantum relaxation due to strong chaos is reported by Prosen p3 
in a perturbative evaluations of fidelity, which is an overlapping integral between the two 
states evolved by slightly different Hamiltonians. As is discussed in p4[, it is hopeful that 



the suppression of quantum relaxations due to strongly chaotic dynamics will have various 
application. In particular, our scenario, which suggests an approach of the dynamical inhibi- 
tion of entanglement, will also provide important applications to quantum communications 



and computations, which require the protections against decoherence M, 24 
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